RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D ; No. 8, August 1965 

Earth-Ionosphere Cavity Resonances and 
the Propagation of ELF Radio Waves 1 

James R. Wait 
National Bureau of Standards, Boulder, Colo. 

(Received February 3, 1965) 

This is a tutorial exposition of the theory of ELF propagation in the earth-ionosphere waveguide. 
Derivations are not given; instead, the plausibility of the results and the physical aspects of certain 
special cases are discussed. A simplified treatment of the Schumann cavity resonance phenomenon 
is also given. 

1. Introduction 

This is intended to be an introduction to electromagnetic wave propagation in the frequency 
range below 3000 c/s and above about 1 c/s. This range has become known as extremely-low- 
frequency (ELF) to distinguish it from the very-low-frequency (VLF) band which is now regarded 
universally as the range 3 kc/s to 30 kc/s. Attention is confined to propagation in the earth-iono- 
sphere waveguide as this is by far the most important mechanism for the transfer of energy at very 
great distances from the source. Included in the discussion is a simplified treatment of the earth- 
ionosphere cavity resonance phenomenon with special reference to the pioneering work of W. O. 
Schumann. The influence of nonuniformities and surface irregularities on the cavity resonances 
is also considered in an appendix. 

References to past work in this field will be made at the appropriate place in the paper. How- 
ever, at this point, the existence of two earlier review articles is mentioned. One is a chapter in 
a book [Wait, 1962], while the other is a review paper given at a NATO conference over a year 
ago [Galejs, 1963]. In both these articles, extensive references to related work are given. In 
addition, there has been a very extensive bibliography, prepared under the auspices of NATO, 
by Brock-Nannestad [1962]. Therefore, in the present review, only material is referenced which is 
actually used in the presentation. 



2. Formulation oi the Waveguide Problem 

In this section, the salient features of the conventional theoretical approach are given. The 
model is a homogeneous conducting spherical earth of radius a surrounded by a concentric iono- 
spheric shell of inner radius a + h. Thus, h is the height of the ionosphere above the earth's 
surface. In what follows, it is convenient to introduce the usual spherical coordinate system. 
Also, unless stated to the contrary, the fields are assumed to vary as exp (icot). 

Using the above model, and assuming a vertical electric dipole source located at = 0, r=a, 
the expression for the vertical electric field, at r=a, is 

Idsrj " KjH-1) p t 

E r = . o;f , 2 y 8 n — Pv{— COS 0)' (1) 

ilkha 1 f* Q sin pit 



Paper was presented at the UI.F Symposium, Boulder, Colo., August 17-20, 1964. 
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Here, k = a)/c is the free-space wave number; 8 n is an excitation factor discussed below; / is the 
average current in the source dipole; ds is the length of the source dipole; -n, the intrinsic impedance 
of the air space, = 12GV; P v {~ cos 6) is the Legendre function of argument — cos 6 and (complex) 
order v. 

The complex quantity v contains the crux of the propagation problem. It is related to the 
propagation factor S n for rc = 0, 1, 2, ... by the relation 



v+ 1/2 = kaSn- 



(2) 



The factor S„ is determined from the boundary conditions of the earth ionosphere waveguide. 
At ELF, numerical values of S n may be determined on the basis of a flat-earth model. Assuming 
for the moment that the earth and the ionosphere are homogeneous isotropic media, it follows 
from previous work [Wait, 1962] that 



s„= 1 



irn 
~kh 



2 1 

4 






y, 



(3) 



/here 



NiN g ' 



(4) 



in terms of the refractive indices, N< f and Nu of the earth and the ionosphere, respectively. This 
equation is valid subject to the condition |A|M < < 1. The sign of the radical is chosen such that 
the real part of S n is positive and the imaginary part of S u is negative. 
For the important case n = 0, (3) simplifies to 



So- 



1 



l kh 



1/2 



(5a) 



Furthermore, since | A | khl(7rn) 2 « 1, when n > 0, it follows that 

1/2 



S n = 



1- 



kh 



2 .2A 

— 1 



kh 



(5b) 



for n=l, 2, 3, ... . 

When the ionosphere is a homogeneous isotropic medium of conductivity en and the earth 
is a homogeneous isotropic medium of conductivity a g , it is well known that 



Thus 



or, in most instances at ELF, 



N\ = a-i/tieou) and N* = o- g /(ie (o). 



wffl'fa& 



—«©'"• 



because cr g >> a,. 
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It might be mentioned in passing that this simple expression for A is only valid for the idealized 
model discussed here. However, it is interesting to note that a horizontally stratified ionosphere 
will give rise to a representation of the form 



A-(l + fl(ifc?, 



where the complex function Q has been tabulated and graphed for a wide variety of models [Wait, 
1962]. 

3. Peripheral Waves 

The physical significance of (1) may now be discussed. The Legendre function is well ap- 
proximated by 



/y-cos 6) = ( ^— V 'cos [(i> + ^)(7r-0)--(7r/4) 

\ttv sin 6/ |_\ 2/ 



(6) 



if | v\ » 1 and is not near or rr. Thus, the individual terms (i.e., modes) in the series in (1) 
are simply proportional to 

, . * 1/2 cos [kaSn(7T -0)- (tt/4)] 
(sin 0y lz 

which, apart from a constant factor may be identified as the linear combination of two peripheral 
waves of the form 

(sin eyi* exp k~ ikaS * e 1> 

and 

— — — jTj exp [— ikaSn(27T — 0)] exp (i7r/2). 
(sin By 11 

It is apparent that these are waves traveling in opposing directions along the two respective great 
circle paths of respective length ad and a(27r — 6). An interesting feature is the 7r/2 phase advance 
which the wave traveling on the long path accumulates as it passes through the antipodal point 
6 = tt. This, of course, is a characteristic feature of symmetrical waves passing through a perfect 
focus. It is seen that the phase velocity of these waves or modes of order n is c/ReS,, where c 
is the velocity of light. Because ReS n > 1, the waves are "slow" relative to unimpeded plane 
waves in free space. The attenuation is given by — klmS n in nepers per unit distance apart from 
geometrical effects which are contained implicitly in the multiplicative factor l/(sin 0) 1/2 . 

The linear combination of the two peripheral waves forms a standing wave pattern whose 
distance 6 W , between minima, is approximately given by 

kdtnReSn = rr or d,„ — A./ (2 Re S n ) 
where A. is the free-space wavelength. 

4. Antipodal Effect 

The asymptotic approximation given above for P v (— cos 6) is not usable when 6 is very near rr. 
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However, in this region, a modified form shows that the field is proportional to the Bessel function 

Jo[kaS n (7T-0)l 
which is finite at 6 = it yet still behaves as 

(sin O)- 1 ' 2 cos [kaS n (TT -6)- (tt/4)] 
when the argument is large compared with unity. 

5. Near Field Behavior 

At very short distances (i.e., 6 near 0), there is another expansion to the Legendre function 
which is quite useful. In the present notation, this result is given by 

PA ^™* 6) s Hf\kS n p) - ± g) 2 H^(kS n p)+ terms in (pla)\ (p/af, etc., (7) 

where v+ h = kpS n and p = a0. When the great circle distance p is reasonably small compared 
with the earth's radius a, only the first term is important. 

6. General Field Representations 

For convenience in what follows, it is desirable to express the field components as a ratio to 
the quantity 

E = i(r)l\)Ids(e- ik e)lp, (8) 

which is the radiation field of the source at a distance p on a flat perfectly conducting earth. With 
this normalized factor, it is not difficult to show that 

E r = E Q W, (9) 

where 

W Wtiplh)** J A n SlH^(kSnP), (10) 

where terms containing (p/a) 2 , (p/a) 4 , etc., have been omitted. An explicit formula for the excita- 
tion factor A n is given by [Wait, 1961], 

AOn^n 
n = ~ 



C 2 -A 2 



where 8 n is obtained from 



In this equation, 



On 



dmoRjovdc 



2khRi(C)R g (Q J 



(11) 



^(C)- c+v **- Ng > 
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and 

The complex factors C n are solutions of the modal equation, which for ELF may be well repre- 
sented by its flat-earth version. Explicitly, it is given by [Wait, 1962] 

R g (C)Ri(C) exp (— ilkhC) = exp (— iTtin), Q2) 

for /i = 0, 1, 2, ... . An approximated solution is given by (3) above or the simplified forms given 
by (5a) and (5b). 

7. Properties of the Modes 

The behavior of the modes in the earth-ionosphere waveguide is best understood by an exam- 
ination of the complex function W defined by (10). When kp > > 1, corresponding to the "far- 
zone," it follows from (10) that 

(p/X) 112 
W ~ hl\ 6XP [ii27T P /k) " lTr ^ 2 A » S n /2 exp [- i2irSnplk], (13) 

n=() 

where A is the free-space wavelength. Since E r = E W, it follows that the vertical electric field 
of each of the modes is proportional to 

(Up) 1 ' 2 exp [(2irpl\)ImS n ] exp [- i(27rp/X)ReS»] , 

which have the character of outgoing cylindrical waves with damping according to — (27r/k)lmS n 
nepers per unit distance. Actually, at ELF in the far zone, only the n = term is really important 
because of the evanescent character of the modes of order one and higher. 

At short distances, where 2rrpl\ is comparable with unity, some of the evanescent modes are 
important. For this region, it appears that the only recourse is to sum the Hankel function series 
given by (10). Some results, using this approach, have been reported previously [Wait, 1962; 
Wait and Carter, I960]. At very short distances, where p << h, it becomes possible to neglect 
the presence of the ionosphere altogether. In this case, 

w ^( l -t P -ik)^~ik» (14) 

which shows that the vertical electric field varies according to 1/p 3 when (IrrpiX) is small compared 
with unity. 

8. Properties of the Excitation Factor 

Some comments regarding the excitation factor are in order. By carrying out the differen- 
tiations indicated in (11), it follows that 



A n = 



Cl-M 



l | lAy | ih 



kh(C*-Al) *A(CJ-A?)J 



(15) 



which is quite general provided that A* and A^ are independent of C. In addition, it is assumed 
that |A 7 | < < |Aj| which is certainly valid at ELF. Equation (15) then simplifies to 
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An = 



1 . sin 2khCn (16) 

IhhCn 



where use has been made of the mode equation (12). A further simplification, which is sometimes 
used at ELF, corresponds to using 

A = 2 since khC < < 1 



A n = 1 for n — 1, 2, 3, . . . since 2khC n = 2mr. 

The latter approximation is not very good unless /?,• has a magnitude near one. However, the error 
introduced by its use is negligible since these modes are evanescent in any case. 

9. Height-Gain Effects 

In the above discussion, the transmitting and the receiving antenna have both been assumed 
to be located on the earth's surface. The' effect of raising either or both of these is readily ac- 
counted for by the use of a height-gain function F n (z). For example, the vertical electric field of 
mode n at height z, is F n (z) times the vertical electric field at height zero. By following through the 
formal derivation or by direct reasoning, it is found that 

exp(ikC n z) -f R g (Cn)exp(— ikC n z) 

Fn(z)= t+rJc;) w 

This obviously has the proper z-dependence, being the superposition of a downgoing and an 
upgoing wave. Furthermore, it has the required property that F w (0) = 1. In fact, a series develop- 
ment shows that 

(kC z) 2 
F n (z) = 1 + iAgkz- - " ) -h terms in (kzf, (kz)\ etc. 



10. Generalized Representation 

A representation for E z which is applicable at very short, intermediate, and great distances 
may be obtained by generalizing (10), taking into account the (sin d)' 1 ' 2 dependence which is also 
consistent with (7). Also, the appropriate height-gain functions are inserted to account for the 
finite height z of the transmitting antenna and the finite height z of the receiving antenna (assuming 
that both are equivalent to vertical electric dipoles). Thus, 

/ V 12 
where 



W iTT( P lh)e ik " J \ n SlH«\kSnp)F n (Zo)F n {z), (19) 

M = 

where p — a6. This result should be valid for all angular distances 6 up until 6 is near tt. The 
heights z and z are to be restricted to the range to h. 

For the same idealized model as described above, the horizontal electric field E$ of the vertical 
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dipole source may be written 



Eq — — E{)S 



o V' 2 



sin 



yhere 



Trp 



S= h eiKP 2 AnSnH ^ (kSnp)Fn(Zo)Gn(zl 



(20) 



(21) 



with 



(^n\Z) — C n 



exp jikCnz) — R(j{C n ) exp ( — ikC n z) 
l+R g (Cn) 



It is also found that the magnetic field, which has only a (f> component, is given by 

6 v/ 2 



Hj, — — E()T 



sin 



(22) 



where 



TTp 



T = ~h e ' k " 2 ^nS n m\kS nP ) F«(z )F n {z). 



(23) 



In summary, we see that the far-zone expressions for the field functions W , S, and T may be 
written in the compact form 



W 

S 



.(p/X) 



1/2 



(A/A) 



g«jrp/X e -fcr/4V A B e _f2,rS n i 



p/A 



S;fF„(2) 
-S,fC„(z) 

S!/ 2 f„( 2 ) 



F„(z ), 



(24) 



which is valid for (p/X)> > 1. 

It is important to note that, quite generally, the ratio S to T, for z = 0, is A fl . This is a necessary 

consequence of the boundary condition Ee = —r)^ g H ( t ) on the earth's surface. It is also interesting 
to note the ratio W/T or —E r l(r)H<t>) is only approximately unity, being equal to S n in the far zone. 
At the shorter distances, when p < X, the ratio W\T will depart significantly from unity. Calcula- 
tions which illustrate this phenomenon are available in the forms of curves of W\T as a function of 
distance [Wait and Carter, I960]. 

It is possible to excite ELF electromagnetic waves by other than a vertical electric dipole. 
For example, a horizontal electric dipole of moment Ids and oriented in the (/> = direction will give 
rise to a vertical electric field E r according to the relation [Wait, 1962] 



6 Y /2 
Er = E S[— — -I cos </>, 
x sin 6/ 



where 



TTD 



5 A^ P 2 knSnH?KkSnp)Gn(zo)F n (zY 



(25) 



(26) 



The function S has the same form as S defined by (21) except that z and zq are interchanged. 
The function G n (zo) is defined by (22) which is valid for Zo in the range to h. Of some interest is 
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the case where the horizontal electric dipole is buried to depth z = — z\. Then, for the same dipole 
current moment, G n (zo) is replaced by G n {zi) where 

Gnizi) = exp [- i(27T Zl lk) (N* - S^ 12 ] = exp [- i(2irzil\)N g l (27) 



where N g = ( — y=r) (— -J is the refractive index of the earth which is assumed homogeneous. 

The horizontal electric dipole also radiates horizontally polarized waves even when the iono- 
sphere is isotropic. As indicated previously [Wait, 1961], this gives rise to a radiated H r com- 
ponent given by 



r)H r = -sin<l>EoTn, (28) 



where 



where 



irp 



T « = h eikP 2 AfM 2 W mP )FtW/M, (29) 



C 2 1 

l ^ kh[Cl-{^fY kh[Cl-{Mlf] 

in analogy to (15) for A m . In the case of the horizontally polarized waves, the CVs are solutions 
of the modal equation for TE (transverse electric) modes. It is given by 

R h g (C)R^{C) exp ( - ilkhC) = exp ( - 2mn\ (31) 

where 

9 C+A* d l C + A? 
The appropriate height-gain function is given by 

FJZ) T+R^) (32) 

for < z < h. 

While the attenuation rate for the horizontally polarized modes, above cutoff, is quite low, 
the excitation factor A* is very small as indicated by the approximation 

Ai^-ClK^f^-Ald. (33) 

On the other hand, the height-gain function is a rapidly increasing function of height. Thus, 

Fife) = 1 + ikz/A g - (kC m z) 2 /2 + . . . , (34) 

where, as usual, k=2irlk. Thus, for a highly elevated horizontal electric dipole, the received 
horizontally polarized field may be appreciable. 
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1 1 . Cavity Resonances 

Within the framework of the present theory it is interesting to see how the Schumann reso- 
nance phenomenon emerges. The starting point is the expansion formula [Schumann, 1957]. 

PA-x) = 1 - 2/i+l 

sin pit tt^ nW rtn+l)-v(v+l) (35) 

where the summation is over integral values of n. Introducing this result into (1) readily leads to 
the result 

„ ildsv(v+l) 2y n , , 2/i+l 

Er= A \ ' y Prix) -—_ 7TT — ~ T7 (36) 

QTTa^ecoh *^q n{n-r 1) — v{v+ 1) v ' 

where x = cos 0. The early terms of the series are then proportional to 

Po(*) = l, 

Pi(#) = cos fl, 

P 2 (*)=s(3cos 2 0-1). (37) 

The 71 = term corresponds to quasi-static field configuration and it is of no particular interest in 
the present discussion. The terms corresponding to n=l, 2, 3, . . . are directly related to the 
Schumann cavity resonances. This may be demonstrated in a fairly simple manner. For a 
homogeneous isotropic ionosphere, 

1 /feftA 1 ' 1 

A '- = /vTUr) (38) 

in terms of the conductivity or, of the ionized medium. Then, using the relation of the operational 
calculus, ico is replaced by p, which leads readily to the relation 

(c/a) 2 K^+l) = ~p 2 -p 3/2 a, (39) 

where a = l/[/i(cr;/x) 1/2 ]. Assuming that the current moment is an impulsive function of time, 

Ids = (Ids) d(t), (40) 

where 8(0 is the unit impulse function at £ = 0. The Laplace transform of the source moment is 
then given by 



Idse-t )t dt = (Ids)o' (41) 

It is now a simple matter to show that the transform of the field is given by 

(Ids) . * 2/i+l 



a>l=(cla) 2 n(n+ 1), 
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(42b) 



and 



~ 4>7ra 2 € 

Ce = — — (42c) 



The actual time response e r (t) is now found from 

EAp)= \%At)e-*dt 9 (43) 

which may be inverted using standard techniques in operational calculus. 
For small damping, it is found [Wait, 1962] that 

e r (t) = { j^ J Pn(x)(2n + l)e-<V C os a#, (44) 



and 



w,; = w w (i- 2 3 /2wl/2 j, 



OUli ' 



2 3/2' 

This transient result is valid for a 2 \ > > 1 or t > > 1/(A 2 ctj/x) and, therefore, it is only usable for suffi- 
ciently small damping. Actually, in the case of zero damping, a = 0, and o)n — 0) n , which includes 
the static field for a>o = 0, and the cavity modes o>i, 0)2, o>3, .... For a = 6400 km, 

o)i/2tt=10.6 c/s, 

C0 2 /277=18.3 c/s, 

o) 3 /2t7 = 25.9 c/s. 

The finite conductivity of the ionospheric shell tends to reduce these frequencies by about 20 
percent. 

12. Final Remarks 

In the present theoretical exposition, a very simple model has been assumed, namely, a 
spherical homogeneous earth with a concentric homogeneous isotropic ionosphere. In recent 
years, there has been much effort devoted to more complicated models. Since this is a short 
review paper it is not possible to develop these generalizations and extensions. However, for the 
sake of completeness, the required modifications to the simple theory are briefly mentioned. 

The presence of the earth's magnetic field renders the ionosphere anisotropic. Fortunately, 
when the earth's magnetic field is steeply dipping it is appropriate to invoke the QL (quasi-longi- 
tudinal) approximation. With this idealization, it is possible to say [Wait, 1962] that for an 
individual mode, 

Attenuation rate with magnetic field _ cos (r/2) 
Attenuation rate without magnetic field (cos r) 1/2 

. (Ol Longitudinal component of gyrofrequency 

where tan r = — = ^ 1V . 7 • 

v Collision frequency 

While this appears to be an exceptionally simple prescription, it does tally with current ideas on 
the subject as discussed by Galejs [1963] in his excellent article. It is seen that for r > 0, the 
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effect of the earth's magnetic field is to increase the attenuation and, consequently, the damping 
of the Schumann resonances are increased. 

Finally, one should mention again that the nonsharpness of the lower edge of the ionosphere 
will play an important role in ELF propagation [Wait, 1962]. As discussed by Galejs [1963], 
this will modify considerably the frequency characteristic of predictions based on the simple 
homogeneous model. Another factor considered by both Fligel [1962] and Galejs [1963] is the 
influence of heavy ions which will modify the effective refractive index of the medium to some 
extent. 



13. Appendix. Application of Perturbation Theory to the Nonuniform Earth- 
Ionosphere Cavity 

In the main text of this paper, and in all published theoretical work on the Schumann resonance 
theory, the earth and the ionosphere are assumed to be concentrically stratified. In this appendix, 
it is proposed to consider the effect of nonuniformities of the ionosphere. In particular, by adapt- 
ing a rather simple model it is shown how the resonant frequencies of the cavity could be modified 
by localized nonuniformities. The theory is somewhat analogous to that developed for microwave 
cavities under the influence of wall imperfections [Harrington, 1961]. 

We consider a concentric spherical cavity of inner radius a and outer radius a + h. The 
surface impedance of the inner surface is zero, while that of the outer surface is Z, which is taken 
to be constant everywhere. The resonance frequencies for this situation are denoted o) n where 
n=l, 2, 3, . . . and the associated electric and magnetic fields within the cavity are E and H, 
respectively. We now consider that the surface impedance Z of the outer surface is changed to 
Z' such that the difference Z' —Z may be a function of and (/) in the spherical coordinate system. 
The resonance frequencies are now changed to (o' n . The associated electric and magnetic fields 
are E' and H', respectively. 

Maxwell equations for the unmodified cavity are: 

-VXE = i/JLCOnH, (Al) 

and 

VxH = i6w„E, (A2) 

while, for the modified system, 

-VxE'^iKH', (A3) 

and 

VxH' = i€ W ;,E'. (A4) 

By making use of the vector identity, 

V-(AxB) = B-VxA-A-VxB, (A5) 

it is now a simple matter to show that 

V- [H'xE-HxE']=;K-co„)[eE-E'-^H-H']. (A6) 

By integrating (A6) over the volume V of the cavity and making use of the divergence theorem, 
it follows that 



f f [H'xE-HxE] 'iida = i(co' n -o)n) [ [ f[eE-E'-/xH-H'] dv, 



(A7) 
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where n is the outward unit normal vector to the surface A which bounds the volume V. Now, 
the boundary conditions on the outer surface are 

nXE = ZH,, (A8) 

and 

nXE'=Z'H,\ (A9) 

where H, and H,' are the tangential magnetic fields for the unmodified and the modified system. 
Because of the assumed perfect conductivity of the inner surface, that part of the surface integral 
in (A7) over A vanishes. Therefore, by using (A8) and (A9), the general relation (A7) simplifies to 



[[ (Z'-Z)H r ll' t da = i(a>' n -(o) f J J [eE • E'-ptH • H'] dv, 



(A10) 



where the area A is the outer surface of the concentric resonator. 

Equation (A10) is an exact formulation of the problem. Thus, if the solution for the uniform 
cavity is known, this equation is a connecting relation between the fields for the modified cavity 
in terms of the change Z' — Z of the surface impedance of its walls. Unfortunately, it is not an 
explicit formula. Instead, it is a rather complicated integral equation. Nevertheless, if the modi- 
fications in the cavity are small it can be expected that the fields E', H' do not differ appreciably 
from E, H. Thus, the first-order perturbation formula for the change of resonant frequency of 
mode n is given by 

i J J (Z'-Z)Ut-K t da 

ft>; " C °" = J/Jt/iH-H-cE-EJdt; (AU) 

V 

where Z' — Z is the specified change of surface impedance and the field quantities are all assumed 
known. Implicit in the statement of (All) is the neglect of coupling between modes of different 
order. For example, if the uniform cavity is oscillating in a single mode, it is possible for additional 
modes to be excited when the wall impedance is modified in a localized region. 

Here, we shall employ (All) with some rather idealized assumptions about the nature of the 
fields in the uniform or unperturbed cavity. Specifically, we assumed that the unperturbed cavity 
is excited by a symmetrical source (e.g., radially oriented electric dipole) at = 0. From section 
11, it is seen that the fields associated with this mode are then given by 

E r = An((On)Pn(cOS 0), (A12) 

and 

r)H* = - L4 n (ft>»)Pi(cos 0) (—\ (A13) 

\(o n a/ 

where A n (a> n ) is the amplitude of the mode at angular frequency cu n . The other field components 
are either zero or vanishingly small. 

Using (A12) and (A13) in conjunction with (A10), it is found, without difficulty, that 



r \ h 2 C 2n f 77 

— )—A [Pn(cos 0)] 2 sin 6 (Z'-Z)dOd(f) 

n aJ if Jo Jo 

= t€«-a>n) T ( 2U ("I [P»(cos 0)] 2 -h [Pi(cos 0)] 2 ("^Vl r 2 sin OdOd^dr. (A14) 
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Without introducing appreciable error, r on the right-hand side of (A14) may be replaced by b. 
Then, by making use of the identities 



/: 



[P n (cos6)] 2 sin Odd = 



2ra+l 



(A15) 



and 



/: 



[Pi(cos0)] 2 sin 6d6 = 



2tt(tt+l) 
2n + l 



(A16) 



it follows, after some slight rearranging, that 



Utt/J 



_2; 

2n( 



n + 1 1 



fifiV _J_ +1 1 

C / 71(71+1) J 



T27T Ttt 

J J o 8(fl,<^)[Pi(cose)]«(8ine)rfed0, ( ai7) 



where 8(0, $) = (Z' -Z)li). 

A very simple illustration of the perturbation formula is to let Z = 0, which means that the 
unperturbed cavity has perfectly conducting walls. For this special case, we also take Z' inde- 
pendent of and (/>. In other words, 



where A is a constant. Then, by making use of (42b) and (A16), it follows that 

_ *Ac 



(A18) 



where 



(o n — (c/a) [n{n + l)] 1 



(A 19) 



This formula for a)' n is a remarkably simple result. Its validity may be tested by comparing it with 
the result based on a direct solution of the modal equation for the uniform concentric resonator. 
From such a previous analysis, it was shown that [Wait, 1964] 



0) n = G) n 



l-i 



cA 



-1/2 



(A20) 



where (o n is defined by (A19) and A is the (normalized) surface impedance at frequency (d„. It is 
easy to verify that (A20) reduces to (A18) when cA/(co n h) < 1. As expected, the first-order per- 
turbation theory is only valid when the frequency w^, for the modified system, is not too different 
from the frequency o> n of the unmodified system. 

Equation (A17) may be used to study the influence of a small localized disturbance. For 
example, we imagine that Z' is different from Z only over a small region bounded in by O — (80o/2) 
and 0o + (80o/2), and bounded in by $o — (8</>o/2) and c/>o + (8</>o/2). The resonant frequency 
o) n for the unperturbed cavity is to be found from [Wait, 1964] 



: (*>0n 



1 



cA Y" 2 

a>nh\ 



(A21) 



1069 



where coo n = {c/a) [n{n + 1)] 1/2 , and where k = Zjy) is a constant. Solving (A21) as a quadratic in 
Q) n , it is seen that 



(O n — (Oqh 



1- 



1 /cA 



4a4 \ A 



^ 2 icA 
2A' 



(A22) 



which approximates to 






(A23) 



if 



cA 



<t 1. In terms of co„ which is known, the resonant frequency co^ for the perturbed resonant 



2h(t)o n 

frequency co' n is now given by 



M l 


Ct) n — 




:) 


" 2n+l 
2n(n+l)_ 


w n 




i+4 

L WoJ 





8(80 o )(80o)[Pi(cos 0o)] 2 sin 0o, 



(A24) 



where 8 is the average value of 8(0, <t>) over the small perturbed region. The explicit do dependence 
of the right-hand side of (A24), for the first two modes, is evident from the identities 

[P}(cos o )] 2 = sin 2 0o, 

[PJ(cos do) f=9 sin 2 o cos 2 o . 

These show that the first mdde is most affected by surface impedance perturbations near the 
equatorial plane whereas the second mode is most affected by perturbations near O = — 45°. 

It would appear that the perturbation theory of cavity resonators should be a useful approach 
to investigating the influence of localized nonuniformities of the ionosphere on the Schumann-type 
resonances. 
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